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метода конечных элементов и комплекснозначных искусственных нейронных сетей* 
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Рассматривается применение комплекснозначных искусственных нейронных сетей (КИНС) в обратной задаче 
идентификации упругих и диссипативных свойств деформируемого твёрдого тела. Дополнительной инфор- 
мацией для решения обратной задачи являются компоненты вектора смещений, измеренные в наборе точек 
на границе тела (позиционное измерение), совершающего гармонические колебания в области первой резо- 
нансной частоты. Процесс измерения смещений в работе моделируется расчётом в конечноэлементном паке- 
те АМ$У$, построением амплитудно-частотных характеристик (АЧХ) смещений и выбором их значений для 
некоторого набора частот (частотное измерение). В приведённом численном примере исследуются вопросы 
точности идентификации модуля упругости и добротности материала в зависимости от числа точек измере- 
ния и их расположения, а также от архитектуры нейронной сети и длительности процесса её обучения, ко- 
торый осуществляется с помощью алгоритма комплекснозначного обратного распространения ошибки (КОР). 
Ключевые слова: комплекснозначные искусственные нейронные сети, идентификация механических 
свойств, метод конечных элементов. 


Введение. В настоящее время искусственные нейронные сети (ИНС) [1] широко применяются в 
различных областях науки. В механике они используются для решения коэффициентных и гео- 
метрических обратных задач. Так, применению ИНС в задачах неразрушающего контроля посвя- 
щены работы [2-6], а определению механических свойств твёрдых тел — работы [7-11]. В начале 
90-х годов прошлого столетия в работах Т. Ма [12, 13] были предложены КИНС, которые в 
настоящее время широко используются для решения прикладных задач [14—16]. КИНС, парамет- 
ры которых (веса, пороговые значения, входы и выходы) являются комплексными числами, при- 
меняются в различных областях современной техники, таких как оптоэлектроника, воспроизведе- 
ние изображений, синтез речи, машинное зрение, дистанционный сбор данных, квантовые аппа- 
раты, пространственно-временной анализ физиологических нейронных аппаратов и систем. При- 
менение КИНС в задачах механики является новой областью исследований, которая начала раз- 
виваться только в последние годы. 

В настоящей работе КИНС применяется к решению обратной коэффициентной задачи 
идентификации упругих (модуль Юнга) и диссипативных (добротность) свойств деформируемого 
твёрдого тела. Дополнительная информация для решения этой обратной задачи связана с изме- 
рением поля смещения на границе тела (в дискретном наборе точек), совершающего гармониче- 
ские колебания в окрестности первой резонансной частоты. В приведённом ниже примере иссле- 
дуются вопросы точности идентификации механических свойств материала в зависимости от чис- 
ла точек измерения и их расположения, а также от архитектуры нейронной сети и длительности 
(количества эпох) процесса её обучения, который осуществляется с помощью алгоритма КОР. 
Постановка обратной задачи идентификации механических свойств. Рассматриваются 
установившиеся гармонические колебания прямоугольной области (ах) с круговой частотой 
в рамках плоской задачи теории упругости (рис. 1). Прямоугольник закреплён по левой боковой 
стороне $, , в правом нижнем углу действует вертикальная сила амплитуды А, ‚ остальная грани- 
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ца $, свободна от механических напряжений. Механические свойства материала описываются 
модулем Юнга Е, коэффициентом Пуассона \ и добротностью О. 








#(2;2 
Бе 


Рис. 1. Модель тела с номерами точек измерения АЧХ смещений 


Уравнения колебаний тела имеют вид [18]: 


о; =-р (6 — юа)и,, К, =1,2, 


ю, 
- - 1 
би = Си (1+ В)Е у, 0, = С. (1 + ЮВ) Е1› , к = 5 (м. +И, к). (1) 
Граничные условия предполагают наличие сосредоточенной силы в точке 2 (рис.1) и од- 
нородные условия на остальной границе: 

и, |5, =0, суп, =0. (2) 

Дополнительной информацией для решения обратной задачи являются смещения и, , из- 
меренные в точках 2—5 (рис. 1): 

и =И, +0 т=2,3,4,5, (3) 


где о,,, &, — компоненты тензоров напряжений и деформаций соответственно; р — плотность. 


ти М И + И 
Коэффициенты упругости соответствуют изотропному телу с., =с,, =Л+2и, с, =с,, =Л, 
са =, Л, н — коэффициенты Ламе. Частота колебаний а совпадает с первой собственной ча- 


стотой резонанса для тела, в котором не учитывается диссипация механической энергии. Коэф- 
фициенты а и В, характеризующие диссипацию, вычисляются по методике, приведённой в [17]. 


Архитектура КИНС. В качестве инструмента решения обратной задачи идентификации механи- 
ческих свойств упругого тела в настоящей работе используется КИНС (на рис. 2). 


Х. =Х:. + Мы И =У, +, 








Хи = Хи к 1 У, — 28 ая 7: 


Входной слой Скрытый слой Выходной слой 


Рис. 2. Архитектура КИНС 


Рассмотрим КИНС, состоящую из 3 слоёв: входной слой, скрытый слой, выходной слой 


(рис. 2). 

Набор весовых коэффициентов (И/, И/ ) показывает связи между данным нейроном и все- 
ми нейронами предыдущего слоя. Причём, каждый весовой коэффициент является комплексно- 
значным. Процесс обучения сети состоит в следующем: 

1. Определяется случайная начальная конфигурация весов с равномерным распределением 
действительной и мнимой части в пределах (—1,0, 1,0). 


85 


Инженерное дело, технологии и технические науки 








2. Входной набор данных (Х, =Х,, +, Х, =Х+Х,,Х,=Х,+Х,), на котором сеть 
должна быть обучена, вводится во входной слой сети, затем вычисляются выходные данные. 

3. При сравнении полученных данных и известных выходных данных (для рассматриваемо- 
го входного набора) вводится вектор ошибки, квадратичная невязка между полученными и из- 
вестными (опытными) данными: 

У =УИ/„Х, =Х + =2, (4) 
т 
1 1 


=. т, (5) 


в -(} то" ве) ва(о. «тик, 


где И/„ — вес связи нейронов с номерами пи 17; Х„ — входной сигнал от нейрона т; Е — сиг- 











_ (6) 


мовидная НКЦИЯ; Е — ошибка ДЛЯ ВЫХОДНОГО СЛОЯ, О’ — активность данного узла; То — 
р п 7 п 


требуемое выходное значение. 
4. Вектор ошибки используется при модифицировании весовых коэффициентов выходного 
слоя, величина ошибки уменьшается при повторной подаче того же набора входных данных: 


И № И + ДИ, м и т Н„АА, ' (7) 
АХ, = = Ве [ 5" | (1-Ве[О, |)ве[0, | + 1т(5" |(1 - т[О, ])т[0, ], (8) 
где ^, — порог нейронов входных слоёв сети с номером /; 6” =Т, -О, — обозначение ошибки 


между фактическим и целевым значениями выходных нейронов; Н» — выходы из скрытых 


нейронов. 

5. На следующем шаге весовые коэффициенты скрытого слоя модифицируются аналогично 
предыдущим. Затем сравниваются выходные сигналы нейронов скрытого слоя и входные сигналы 
нейронов выходного слоя. В результате такого сравнения формируем векторы ошибки для скры- 
того слоя: 


И’, =И/, + АИ, = И, + 18, (9) 
де, = в -ве[н, вен, 15. (е[5'" а -веГо, ео, }ве[и/, + 
15" (1 оо, ито, ит [м] - 
Е | =ш[Н,„ тм, | У (ве [6 |(1-Ве[о, |)Ве[о, ]1т[И/, | - 
-ит[5" (1 ито, уия(о, ]ве[и,]]}, 


где И’, — вес между входными и скрытыми нейронами с номерами [и т; Г, — ВЫХОДЫ ИЗ ВХОД- 


(10) 


ных нейронов; 9, — пороги нейронов скрытых слоёв сети. 


6. Сеть обучается путём предъявления каждого входного набора данных и последующего 
обратного распространения ошибки. Такой цикл повторяется много раз. Обучение заканчивается 
при достижении для каждого из обучающих образов значения функции ошибки, не превосходящего 
некоторого заданного значения =, либо после максимально допустимого числа итераций (эпох). 

После достаточного обучения сеть способна моделировать функцию, связывающую значе- 
ния входных и выходных переменных, и такая сеть может быть использована для прогнозирова- 
ния выходных значений, то есть решения обратной задачи идентификации. 
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Применение КИНС к задаче идентификации модуля Юнга и добротности. Для идентифи- 
кации модуля Юнга и добротности используется описанная выше КИНС, в которой эти параметры 
являются выходными данными. Входными данными являются амплитуды смещений (3), измерен- 
ные на поверхности тела. Для обучения КИНС составляется набор входных и выходных данных, 
первым шагом при его формировании является расчёт собственных резонансных частот, затем 
проводится расчёт установившихся колебаний для тела на этих частотах или в их окрестности с 
учётом диссипации механической энергии. 

На рис. 1 показаны точки 2, 3, 4, 5, в которых проводятся измерения амплитуд перемеще- 
ний. В работе натурный эксперимент заменялся численным расчётом в конечно-элементном паке- 
те АМЗУ5. На рис. 3 представлена амплитудно-частотная характеристика вертикального смещения 
в точке 2 в окрестности первого резонанса. 
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Рис. 3. Амплитудно-частотная характеристика смещения (/, , измеренная в точке 2 


В табл. 1 представлены входные и выходные данные для «измерений» в точке 2. 


























Таблица 1 
Пример данных из задач модального и гармонического анализа 
№ входные данные я выходные данные, м 
Е, Па [6] >. [77 \, и; 

1 5х10° 10 206,35 | 482х107 | -5,2х10° | 1,09х10° | -7,3х105 

2 10х10° 20 291,83 | 241х107 | -2,6х10° | 546х107 | -3,6х10° 
З 15х10° 30 357,41 160х107 | -1,7х10° | 3,64х107 | -2,4х105 
4 20х10° 40 412,7 120х107 | -1,3х10° | 2,73х107 | -1,8х105 
5 25х10? 50 461,42 | 964х107 —1х105 2,18х107 | -1,5х10° 
































В случае одной точки позиционного измерения и лточек измерения на АЧХ, КИНС имеет п 
комплекснозначных входных нейронов: Ш, „+\/, и, И 2+, р, и И т+ Ш, п, 
И, п+, п, О, 2+0 О т+№Ю, п 

Составлено Х (200, 300, 500) векторов данных, 90 % из которых используются для обуче- 
ния, 10 % для тестирования. Затем проводятся компьютерные эксперименты, выполненные с по- 
мощью КИНС. Результаты обучения и тестирования представлены в табл. 2-6. 

На рис. 4 представлены результаты тестирования на 50 примерах (идентификации модуля 
Юнга) обученной КИНС с архитектурой «100-20-1» (табл. 6, строка № 1). Кружками изображены 
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данные для тестирования Е, треугольниками и квадратиками изображены прогнозируемые дан- 
ные Е. , полученные с помощью КИНС. Относительные ошибки д, , б, вычислялись по формуле 


в. Е. -ЕИЕ, 6, 10. - СМ. 


























































































































Таблица 2 
Результаты обучения и тестирования КИНС со 2-й точкой позиционного измерения 
Количество | Архитектура Количество 
[9 
№ Данных К КИНС Эпохи | Ошибка | Точность | Примечание ноЧенньАЧХ 
1 200 2—4-—1 2000 0,073 84,80 % В точке 2 1 
2 200 6-3—1 2000 0,066 89,55 % В точке 2 3 
3 200 100-10-1 2000 0,017 96,29 % В точке 2 50 
Таблица 3 
Результаты обучения и тестирования КИНС в различных точках 
(с одинаковым количеством эпох) 
№ ОЛИЧЕСТВО: Архитектура Эпохи | Ошибка | Точность Примечание 
данных КИНС 
1 200 100-10—1 2000 0,017 96,29 % В точке 2 
2 200 100-10—1 2000 0,018 92,86 % В точке 4 
3 200 100-10-1 2000 0,021 93,78 % В точке 5 
Таблица 4 
Результаты обучения и тестирования КИНС 
с различным количеством данных обучения 
№ ОЛИЧЕ СВО, | АрХИТеКтУра Эпохи | Ошибка | Точность Примечание 
данных КИНС 
1 200 100-10-1 2000 0,017 96,29 % В точке 2 
2 300 100-10-1 2000 0,0057 97,64 % В точке 2 
3 500 100-10-1 2000 0,0053 98,37 % В точке 2 
Таблица 5 


Результаты обучения и тестирования нейронных сетей с различным количеством 
скрытых слоёв при позиционном измерении во 2-й точке 



























































№ ОЛИНЕСТЕО у ирхиектира Эпохи | Ошибка | Точность | Примечание времЯ, 
данных КИНС с 
1 500 100-5-1 2000 0,0121 98,13 % В точке 2 96 
2 500 100-10-1 2000 0,0053 98,37 % В точке 2 219 
3 500 100-20-1 2000 0,0019 99,09 % В точке 2 300 
4 500 100-50-1 2000 0,0067 97,71 % В точке 2 834 
Таблица 6 
Результаты обучения и тестирования нейронной сети с различным 
количеством эпох при позиционном измерении во 2-й точке 
№ т а ра Эпохи | Ошибка | Точность | Примечание ( а 
1 500 100-20-1 2000 0,0019 99,09 % В точке 2 300 
2 500 100-20-1 5000 0,0012 98,18 % В точке 2 784 
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Номер образца теста Номер образца теста 


Рис. 4. Результаты тестирования (Е, ©) при пороге ошибки 10 % 


На рис. 5 представлены результаты тестирования на 50 примерах (рис. 5, а) (идентифика- 
ции модуля Юнга и добротности) обученной КИНС с архитектурой «10-20-1» (табл. 6, строка 
№ 1) и зависимость относительной ошибки (6) от числа эпох (рис. 5, 6). 
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а) 6) 
Рис. 5. Результаты тестирования обученной КИНС с архитектурой «10-20-1»: 
а) результаты прогнозирования, полученные в системе; 6) график ошибки обучения 


Выводы. В результате проведённого исследования разработан метод идентификации упругих и 
диссипативных свойств твёрдого деформируемого тела, использующего данные гармонических 
колебаний на резонансной частоте на основе сочетания метода конечных элементов и КИНС. 
Осуществлена программная реализация метода и проведено его тестирование при использовании 
первой резонансной частоты колебаний. В результате численного эксперимента выявлены архи- 
тектуры КИНС, дающие лучший результат идентификации, а именно: 100 (входных нейронов) — 
20 (скрытых нейронов) — 1 (выходной нейрон). 

Проведена оценка временных затрат, связанных с обучением КИНС. Разработанный метод 
и программы могут быть использованы для определения диссипативных свойств на различных 
частотах (не только первой), а также для более сложных свойств упругого тела, например при 
наличии анизотропии. 
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ЕГАЗТТС АМО ОТ$$ТРАТТУЕ МАТЕВТАЕ РКОРЕКТТЕ$ ОЕТЕВМТМАТТОМ 9$1МС 
СОМВТМАТТОМ ОЕ ЕЕМ АМО СОМРЕЕХ АВТТЕТСТАЕ МЕУВАЕ МЕТМОВК$* 


А. М. боомуем, Мдцуеп Бицу Тгиопд С!апд 


Тре арр/саНоп оЁ те сотриех агий<а/ пеига! пе иогк$ (САММ) {0 {Ве туегзе аепИйсавоп ргоМет оЕ {Ве е/азИс 
апа @храйуе ргорегйез о! деюгта Ме ос (5 сопз!аегеа. Тре аа@#опа!/ пЮгтайоп Ю {пе туегзе ргоМет 1/5 
сотропепёз оЁ пе @расетепЕ уесюг теазигей т а 5ЕЁ оЁ ротё5 аЁ те о Боипаагу (розопа!/ теазиге- 
теп®). ТР зо! регопт5 ваптогис озс!айоп5 т пе ИЕ гезопапЕ тедиепсу. ТПе ргосез$ оЁ @р/асетепЕ теа<- 
игетепЕ 15 ятиаей изтд йе сасшайоп оЁ Ппйе еетепё 5оПуаге АМ$У$, {те Бий4та о! {Те атридЕе- 
Недиепсу свагас{ейИс$ (АЕС) о! пе @расетепЕ апа о! {те 5еесйоп оЁ Те! уа/иез Юг а $ЕЁ о! Педиепсие5 
(Педиепсу теазигетеп®). т {Те дмеп питептса! ехатр/е, ргоМетз оп {Те ассигае «епИйсайвоп оЕ {те е/а$Ис 
тоаШиз, апа тейепа! диаМу аерепатд оп {1е питрег оЁ теазиге рот апа {те осайоп, аз ие! аз оп Те 
пеига/ пеймогк агспПесиге апа {те 1епоЁ? ог {пе ватта ргосез$ регогтеа Бу пе сотриех-уа/ие еггог Баск 
ргорадаНоп (СВР) а/дотЁт аге пуезИдаёеа. 

Кеуигога5: сотр/ех-уаие агйса!/ пеига/ певмогк$, Гепийсавоп оЁ теспагиса/ ргорегиез, Ппйе еетегЕ 
тетод. 





* Тне гезеагсй 15 допе ми {пе рагЧа! Япапса! зиррог тот ВЕЕТ (дгагз поз. 13-01-00196-а, 13-01-00943-а). 
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